A hyperons in 2 AGeY Ni + Cu collisions 



EOS Collaboration 

M. Justice^'\ S. Albergo^ F. Bieser^, F.P. Brady^, Z. Caccia^ 

D. A. Cebra^, A.D. ChaconS, J. L. Chance^, Y. Choi^'^, 

S. Costa^ J.B. Elliott^ M.L. Gilkes^ J. A. Hauger^ 

A.S. ffirsch^ E.L. Hjort^ A. Insolia^ D. Keane% 

bo"; J. C. Kintnerb'^ M. A. Lisa^'^ H. Lm% H.S. Matis^, 

g; \ R. McGrath^ M. McMahan^^, C. McParland'^, D. L. Olson^^, 

M.D. Partlan^, N. T. Porile^ R. Potenza^ G. Rai^, 

J. O. Rasmussen^, H. G. Ritter^, J. Romanski'^'^, J. L. Romero^, 

G. V. Russo^ R.P. Scharenberg^ A. Scott^'^ Y. Shao^'^ 

^ ■ B.K. Srivastava^ T. J. M. Symons^, M. TincknelP, C. Tuve^ 

(N ■ 
> 



O 



s. Wang^, P.G. Warren^, D. Weerasundara^'^, H. H. Wieman^, 



g: T. Wienold^'9 K. Wolf^ 

O . 

00 ■ '^Kent State University, Kent, OH 44242 



b 



University of California, Davis, C A 95616 



■ ^ Universitd di Catania and Instituto Nazionale di Fisica Nucleare-Sezione di 

^ . Catania 95129 Catania, Italv 

(D . 

Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 
o; 94720 

^ . ^Purdue University, West Lafayette, IN 47907 

^ State University of New York at Stony Brook, Stony Brook, New York 11794 

X ■ ^Tezas A&M University, College Station, TX 77843 

M'. 

Abstract 

A sample of A's produced 'm. 2 A GeV ^^Ni + '^^^Cu collisions has been obtained 
with the EOS Time Projection Chamber at the Bevalac. Low background in the 
invariant mass distribution allows for the unambiguous demonstration of A directed 
flow. The A nriT spectrum at mid-rapidity has the characteristic shoulder-arm shape 
of particles undergoing radial transverse expansion. A linear dependence of A mul- 
tiplicity on impact parameter is observed, from which a total A + S'^ production 
cross section of 112 ± 24 mb is deduced. Detailed comparisons with the ARC and 
RVUU models are made. 
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1 Introduction 



Relativistic nucleus-nucleus coUisions provide a unique opportunity for study- 
ing hot dense hadronic mattcr. At beam energies below 2 AGeV the focus has 
been on extracting information about the nuclear matter equation of state 
from measurements of the collective flow of nucleons and hght fragments [1]. 
As the beam energy is increased, however, a complete description of the EoS 
should also include strange particle degrees of freedom. Recent theoretical 
studies have indicated that the in-medium properties of strange particles in 
nuclear matter at high densities and temperatures may be constrained by ac- 
curate measurements of their yields, spcctra, and flow in heavy-ion collisions 
[2,3]. Preliminary experimental evidence of a directed flow signal for lambda 
particles has been provided by the EOS coUaboration [4] and by the FOPI 
coUaboration [5]. In this Letter we present the completed analysis of A side- 
ward flow in the EOS data. In addition we present for the first time the EOS 
results on the total yield, the yields as a function of centrality and rapidity, 
and the collective transverse expansion of A's. 



2 Experiment and data analysis 

A fuU description of the EOS experimental setup can be found in Ref.s [6-8]. 
The heart of the setup is a Time Projection Chamber (TPC) which provides 
continuous three dimensional tracking and particle Identification for particles 
with Z < 8. The EOS TPC has a rectangular geometry and operated in anti- 
parallel 1.3 T B and 120 V/cm E ficids for this experiment. The target was 
situated as close to the TPC as possible, approximately 15 cm upstream of the 
first TPC pad row, in order to maximize acceptance. Because this placed it in a 
high magnetic field, copper rather than nickel was used as the target material. 
The ^^Ni beam supplied by the Berkeley Bevalac had an energy of 1.97 A GeV 
at the center of the target. The hardware triggcr consisted of beam defining 
counters located some distancc upstream and a trigger counter located just 
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downstream of thc target. The discriminator threshold on the downstream 
counter was set to veto the ~25% most peripheral events. This trigger cutoff 
produces neghgible bias in the data sample since, as demonstrated below, A's 
are predominantly produced in central events. 

A's are reconstructed through the charged particle decay: A — >■ p + tt", which 
has a branching ratio of approximately 64%. After all TPC tracks in an event 
are found and the overall event vertex has been determined, each pair of pvr" 
tracks is looped over and their point of closest approach is calculated. Pairs 
whose trajectories intersect at a point other than the main vertex are fit with 
a VO hypothesis from which an invariant mass and momentum are extracted. 
The A candidates are then passed through a neural network filter [9] in order to 
reduce the combinatoric background. The resuhing invariant mass distribution 
is shown in Fig. 1. Cutting around the peak in the distribution (1112 MeV/c^ 
< Ma < 1120 MeV/c^) results in a sample of 1797 A candidates, of which 
around 40 are estimated to be background. 

Acceptance correction factors have been calculated by passing a sample of 
7.2 X 10^ minimum bias 2 A GeV Ni + Cu ARC events through a detailed 
GEANT simulation of the EOS detector system and then through the same 
analysis chain as actual data. The ARC model [10] has been successful at 
reproducing a wide range of inclusive observables at AGS energies and it has 
been used to study directed fiow in Au + Au collisions at Bevalac energies [11]. 
The parameterization of the TPC response in the EOS GEANT simulation 
code is similar to the parameterization in the FST simulation program used 
by the STAR coUaboration [12]. Hit merging, spatial resolution, and ADC 
resolution parameters were empirically adjusted to match the data. The over- 
all width and background of the reconstructed Monte Carlo A's are in good 
agreement with the data of Fig. 1. 

Azimuthally integrated A efficiency factors are calculated on a 15 x 25 grid in 
{y, uit) space from the formula: 



/ v Nrec{yi,mT,j) 

e{yi,mT,j)^— — ^, 1 



where A^^ec and Nmc are the reconstructed and input number of A's, respec- 
tively, in the given bin. In principle the efficiency factors also depend on the 
multiplicity or impact parameter. Generating a sufficient number of GEANT 
events to calculate three dimensional efficiencies is impractical in the present 
case, however. Therefore, the approximation has been made that the efficiency 
factorizes into two pieces: e(6, |/,mr) = e{b) ■ e{y,mT)] with the e(6)'s being 
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given by: 



averaged over all values of y and m^- In central and semi-central events track- 
ing efficiencies are independent of impact parameter for the relatively light 
Ni + Cu system (maximum TPC multiplicity of around 70 tracks). In periph- 
eral collisions, howevcr, thc prcsence of a forward going high Z fragment creates 
a hole in the TPC acceptance near beam rapidity. Although the factorization 
assumption is not vahd for such events, the overall systematic error introduced 
is expected to be small, as the hole occurs in a region where there are few A's. 
At mid-rapidity, where the majority of A's are produccd, thc cfhcicncy factors 
are very flat versus pr over the entire range of impact parameters. 

The overall acceptance for A — )• p + tt", averaged over b, y, and pr, is calcu- 
lated to be about 16%. Approximately half of the loss is due to geometrical 
acceptance. Determination of the precise fraction decaying outside the solid 
angle of the detector is facilitated by the fact that, to first order, ARC repro- 
duces the shape of the experimental distribution reasonably well. To second 
order, however, a small systematic uncertainty of < 10% might be introduced 
by the extrapolation of cfRciency factors to regions of phase space where there 
are no counts. Somc A's are also lost through inefficiencies in the track re- 
construcion software. According to the simulations, over ninety percent of the 
A's which decay within the geometrical acceptance are reconstructed. Event 
display comparisons of data with Monte Carlo suggest that the true tracking 
efficiency may bc lower by a few percent. The remaindcr of thc cxpcrimcntal 
inefficiency is due to the cuts used to reducc thc combinatoric background. 
More than 50% of the true A's which are reconstructed get thrown out by 
these cuts. Systematic error in the estimation of the magnitude of this loss 
mechanism is believed to be small. The total systematic error in the final cross 
section from all contributions is estimated to be ten to fifteen percent — half 
the size of the statistical error. 

The open channels for A production in nucleon-nucleon collisions at this beam 

energy arc N N — > KKN (1.58 GeV), N N — > YPKN foUowed by — > A7 
(1.79 GcV), and NN — > KKNti (1.96 GeV). A's originating from E° dccay 
(cr = 2.2 X lO^^cm) are experimentally indistinguishable from primary A's. 
In the ARC event sample approximately 23% of the A's arise from TP decay. 
These A's are included in all of the ARC comparisons shown below. All spectra, 
except for the sideward flow plots of Fig. 5, have also been corrected for 
branching ratios. The dependence of A yield upon centrality produces a natural 
weighting towards central events but no explicit centrality cuts have been 
applied. Attempts to study the effects of centrality cuts on the fiow analyses 
were inconclusive due to insufficient statistics. 
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3 Results 



The cfRcicncy-corrcctcd A center-of-mass rapidity spectrum, normalized to 
beam rapidity, is shown in Fig. 2 along with the ARC comparison. Both spectra 
peak near y = O as expected for this nearly mass-symmetric system. Within 
statistics, the widths of the distributions are similar; howcver, ARC appears 
to overestimate the total yield by approximately 50%. This is not necessarily 
a dcficiency of the cascade assumption, however, as there is a relatively large 
uncertainty in the A cross section in pp colhsions below 2 GeV. The sohd curve 
of Fig. 2 represents the rapidity spectrum of an isotropic thermal distribution 
at the temperature of 106 MeV which is obtained by fitting an exponential to 
the mid-rapidity transverse mass spectrum of the data. The rapidity spectra 
for both the data and the model are considerably broader. This most likely 
indicates an insufficient number of rescatterings for the A's — which have a 
well known forward-backward distribution in elementary pp reactions — to 
become completely isotropic. 

The dependence of A yield upon event centrality is shown in the top panel of 
Fig. 3. The impact parameter scale in this figure was derived by assuming a 
monotonic dependence of total observed charged particle multiplicity in the 
TPC versus impact parameter [13]. The data are well described by a straight 
line fit of (A^a) versus b: (Na) = 0.1361(45)-0.0153(7)6 with axVi^ of 0.6. The 
ARC model, on the other hand, predicts a linear dependence of A yield on the 
number of participants rather than on impact parameter. When plotted versus 
the number of participants, the data bend over and flatten out for the most 
central coUisions. This behavior could be suggestive of a kind of "shadowing" 
effect in which late arriving nucleons from the projectile and target suffer softer 
coUisions with matter already stopped in the coUision zone. Since the initial 
beam energy is only marginally above thrcshold to begin with, A production 
from these late primary nucleons might be suppressed. The fact that the effect 
doesn't also show up in the cascade model, however, tends to cast doubt on 
this explanation. 

The efficiency-corrected mean vr^ multiplicities versus impact parameter are 
shown in the bottom panel of Fig. 3. The 7r~ efficiency is approximately 60% 
and independent of b. The vr^ multiplicities of the data are not as well fit by a 
straight line as are the A's. The (tt") dependence upon participant number is 
reasonably linear, however, consistent with previous observations [14]. Within 
statisitics, the A/tt" ratio is constant as a function of b with a value of 0.010 ± 
0.001. For the most central events, ARC overpredicts the vr" yield by about 
20%. ARC pion multiplicities bend over significantly when plotted versus the 
number of participants. The net result is a strong dependence of the A/tt" 
ratio on centrality in ARC. 
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The straight line fit of thc upper panci of Fig. 3 intersects the abscissa at 
b = 8.9 ± 0.7 fm. Intcgrating the fit out to this vahic of b yiclds a total cross 
section of cta+so = 112 ± 24 mb. As mentioned above, a systematic error of the 
order of 15% from uncertainties in the acceptance corrections should be at- 
tached to this result. For the central 180 mb of the total reaction cross section, 
corresponding to impact parameters b < 2.4 fm, thc cxtractcd cross section is 
Ca+eo = 20 ± 3 mb. This is somewhat higher than, but not in disagreement 
with, the 7.6 ± 2.2 mb cross section for central 1.8 A GeV Ar + KCl collisions 
of Ref. [15] when scaled by the approximately 50% increase in the number of 
participants. The parameterization used in Ref. [3] for the A production cross 
section predicts only a 10% increase in yield in going from a beam energy of 
1.8 A GeV to 2.0 A GeV. 

The mid-rapidity (-0.25 < y < 0.25) A mr spectrum is shown in Fig. 4. In 
this representation, a purely thermal distribution would result in a straight 
line. The data of Fig. 4, however, have the charactcristic shoulder indicative 
of a collective outward expansion depleting the low region [16-21]. Ncarly 
twenty years ago, Siemens and Rasmussen noted that the spectra of protons 
and pions measured at 90° in the center of mass at the Bevalac showed a 
pronounced shoulder at low pr, which they attributed to collective expansion 
[16]. By assuming a spherically symmetric thermalized source at mid-rapidity 
expanding with a constant velocity, f3r, they were able to derive an analytical 
expression which, in terms of transverse mass, is: 



mi drriT 



sinha / T \ T ^ 

7 H cosh a 

V mr/ TTiT 



a V niT/ rriT 



(3) 



where a = Mp/T and 7 = (1 - t^lf^. 

Straightforward application of the Siemens-Rasmussen expression to the data 
of Fig. 4 gives (3r = 0.36 ± .02 and T = 51 ± 1 MeV with a xVj^ of 1.1. 
The interpretation of thc fitted parameters is not so straightforward, however. 
Studies of EOS Au+Au data have shown that differcnt assumptions about the 
expansion velocity profile are equally consistent with the observed spectra. 
In general, the dependence upon the mean radial flow velocity and temperature 
cannot be expressed in simple analytical form as above but must be determined 
from Monte Carlo simulations. In thc casc of the present data no attempt has 
been made to cxtract j3.rS or temperatures for other velocity profiles. The 
fit results should therefore be viewed as providing a quantitative measure of 
collective transverse expansion only within the context of the velocity profile 
assumption underlying Eq. 3. 

Thc rriT spectrum of mid-rapidity protons from the same event sample also 
has a shoulder. An independent fit of the protons results in a /5^ of 0.42 ± .01, 
in reasonable agreement with the /3^ for the A's, with a x^/i^ of 0.75. The 
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proton's "tcmpcrature" , however, is much highcr: 81 ±1 MeV. This difference 
is consistent with the picture of cold A's produccd just above threshold under- 
going an insufficient number of rescatterings to come into thermal equilibrium 
with the hotter protons. 

A two parameter fit of the A's obtained by fixing /3r to zero gives T = 106 ± 5 
MeV with a of 1-9- Fitting the protons with a pure exponential results in 
an inverse slope parameter T = 142 ± 1 MeV with a x^/^ of 1.4. The rapidity 
spectrum obtained from the fit parameters of Eq. 3, assuming a spherically 
symmetric source, is indicated by the dashed curve of Fig. 2. The result is 
practicaUy indistinguishable from a purely thermal distribution at the tem- 
perature of 106 MeV shown as the sohd curve. The ARC A's in Fig. 4 show 
no evidence of a depletion at low ttit and are weU fit by simple exponentials. 
The resulting inverse slope parameter is T = 91 ± 2 MeV. The uit spectra for 
ARC protons also show no evidence for transverse flow. Thermal fits give T 
= 121 ± 1 MeV for fihered and T = 126± 1 MeV for unfiltered ARC protons 
with x^/i^'s of 1.3 and 1.2, respectively. 

Detection of a majority of the charged particles in the TPC, along with the 
presence of directed flow for protons and heavier fragments, allows for the cor- 
relation of A production with the event reaction plane. The Standard trans- 
verse momentum analysis of Danielewicz and Odyniec [22] has been performed 
on the 1797 A events. Details of the analysis along with preliminary results 
obtained from a smaller and less clean event sample are given in Ref. [4]. 
Briefly, a reaction plane for each event is determined from protons and nu- 
clear fragments with Z < 8. The A transverse momentum is then projected 
onto the estimated reaction plane and averaged as a function of rapidity over 
all events. 

The resulting (p^) distributions, normalized to the proton mass, are shown 
in Fig. 5a for A's (filled circles) and protons (open squares). Both have been 
corrected for the 35° dispersion in the estimated reaction plane, as per the 
prescription of Ref. [22]. The conclusion of Ref. [4] has not changed — the A's 
"flow" in the same direction as the protons. The dashed lines of Fig. 5 represent 
straight hne flts to the data over the region: O < y < 0.8. The slopes, (m^/m) x 
d{Px) /dy, ai mid-rapidity are given in the flrst column of Table 1. The fact that 
the A slope is consistent with zero within two Standard deviations is somewhat 
misleading. Backward rapidity points have been excluded from the fits because 
reduced acceptance for < O leads to systematically lower values of \{px)\- 
This effect can be clearly seen in the proton data points, which should be 
nearly antisymmetric about y/vbeam — O from symmetry considerations. The 
efficiency factors used to correct the y and uit spectra have not been applied 
in the case of the event-by-event directed fiow analysis. However, the direction 
of the correction can only be toward larger values of |(pa;)|, which bolsters the 
conclusion of a non-zero positive A flow. 
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The analysis has also been performed on the ARC events. The results. both 
before and after the GEANT fiher, are given in Table 1. Within statistical 
uncertainties the cascade is in agreement with the data for both protons and 
A's. A comparison of columns two and three of the table shows (again with 
large statistical errors) that the experimental acceptance and resolution do 
not distort the extracted flow values. 

Also shown in Table 1 are flow values extracted from the relativistic trans- 
port code, RVUU, by the authors of Ref. [3]. Unhke ARC, where flow effects 
arc produccd by preservation of two-body scattering planes and an unequal 
weighting of attractive and repulsive orbits [11], RVUU incorporates a mean 
fleld potential for lambdas which increases the sidewards flow in the direction 
of the protons. The error bar on the RVUU value in Table 1 spans the range of 
flow values that were obtained in Ref. [3] from a ±20% adjustment in the rela- 
tive strength of the vector part of the mean fleld. Prom the table it is seen that 
the two models produce approximately the same amount of flow for protons 
while ARC A's have, perhaps, slightly more flow. Our flndings are consistent 
with the prediction of Ref. [3] that the magnitude of A directed flow is smaller 
than proton flow. Of course, the hypothesis of A flow and proton flow being 
equal cannot be confldently ruled out with the present statistics. 

The roughly twofold difference between the proton {px/^) of the present work 
and the results of Ref. [5] is dominated by the pT/m > 0.5 selection used in 
Ref. [5]. For the A's, both the transverse and sidewards flow analyses are very 
sensitive to the inclusion of combinatoric backgroTind. The protons from A 
decays have momentum vectors close to their parcnt particle's due to the low 
pion mass. Since non-decay protons, which are dominant in the formation of 
the combinatoric background, have both types of flow, false A's also show flow. 
In order to obtain pure A flow it becomes necessary to either: 1) reduce the 
background to an insigniflcant level as was donc in the present analysis; or 2) 
perform a careful analysis of the flow of the background. Choosing the flrst 
option means sacriflcing some statistics but avoids the inevitable complications 
and uncertainties of method 2. The A invariant mass plot of Ref. [5] includes 
a large combinatoric background. It is unclear whether or how much their 
reported A (p-r/m) are shifted towards their proton (pi/m) by this background. 



4 Conclusions 

In summary, we have obtained a high quality sample of A's produced in 
2 A CcV Ni + Cu collisions. To within 15% the acceptance is well understood, 
allowing efficiency corrections to be made with confldcncc. The A + S" pro- 
duction cross section has a linear dependence upon impact parameter and the 
A/tt" ratio is relatively flat versus impact parameter. A shoulder at low pt 
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in the mid-rapidity A tranverse mass spectrum is most likely due to coUective 
transversc cxpansion. Corrclation of thc A momentum with the reaction plane 
on an event-by-event basis reveals a positive directed flow consistent with the 
predictions of mean field theory. The ARC model reproduces the shape of the 
A rapidity spectrum and generates sidewards fiow in agreement with the data 
but fails to reproduce the shoulders in thc A and proton ttit spectra and the 
dependence of the A yield upon event centrahty. 
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Table 1 

Flow values normahzed to the proton mass in units of MeV/c. 

data unfiltered ARC filtered ARC RVUU 

p 133 ± 10 152 ±4 141 ± 9 140 

A 85 ± 43 163 ± 14 113 ± 31 96 ±17 
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Fig. 1. A invariant mass spectrum. 
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Fig. 3. Mean A and tt multiplicities versus impact parameter. 
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Fig. 5. Average in-plane momentum versus normalized rapidity for (a) data and 
(b) filtered ARC. Open squares represent protons, closed circles represent A's. 
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